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Abstract 

In the context of theories of Kaluza-Klein type, with a large extra dimension, we study self-similar cosmological 
models in 5D that are homogeneous, anisotropic and spatially flat. The "ladder" to go between the physics in 5D 
and 4D is provided by Campbell-Maagard's embedding theorems. We show that the 5-dimensional field equations 
Rab = determine the form of the similarity variable. There are three different possibilities: homothetic, 
conformal and "wave-like" solutions in 5D. We derive the most general homothetic and conformal solutions to 
the 5D field equations. They require the extra dimension to be spacelike, and are given in terms of one arbitrary 
function of the similarity variable and three parameters. The Riemann tensor in 5D is not zero, except in the 
isotropic limit, which corresponds to the case where the parameters are equal to each other. The solutions can 
be used as 5D embeddings for a great variety of 4D homogeneous cosmological models, with and without matter, 
including the Kasner universe. Since the extra dimension is spacelike, the 5D solutions are invariant under the 
exchange of spatial coordinates. Therefore they also embed a family of spatially tnhomogeneous models in AD. 
We show that these models can be interpreted as vacuum solutions in braneworld theory. Our work (I) generalizes 
the 5-D embeddings used for the FLRW models; (II) shows that anisotropic cosmologies are, in general, curved in 
5D, in contrast with FLRW models which can always be embedded in a 5-D Riemann-fiat (Minkowski) manifold; 
(III) reveals that anisotropic cosmologies can be curved and devoid of matter, both in 5-D and 4D, even when the 
metric in 5Z) explicitly depends on the extra coordinate, which is quite different from the isotropic case. 

PACS: 04.50. +h; 04.20.Cv 
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1 Introduction 



In standard general relativity, exact solutions have played a crucial role in the development of many areas of astro- 
physics and cosmology [I] , [2] • Exact solutions (i) provide a route to better and deeper understanding of gravity by 
giving non-perturbative insight into the highly nonlinear gravitational phenomena, (ii) can reveal unforeseen features 
of the theory that might be relevant to more general situations, and (iii) can be used for checking computer codes, 
which is important for the advent of numerical relativity [3] . 

In recent years there has been an increased interest in theories that envision our world as embedded in a universe 
with more than four large extra dimensions. Besides the longstanding theoretical motivations - as, e.g., resolving 
the differences between gravity and quantum field theory, and unifying all forces of nature - the interest in these 
theories has now a practical stimulus. Namely, that any solution of the N-dimensional Einstein equations with source 
can be locally embedded in a {N + l)-dimensional Ricci-flat manifold, which is guaranteed by Campbell-Maagard's 
embbedings theorems [l]-[2]. In particular, this means that any solution of the AD Einstein equations with matter 
G^i/ = SttT^j/ (throughout the paper we set c = G = 1) can always be locally embedded in a solution of the vacuum 
Einstein field equations Rab = in 5D. 

Finding an embedding for a particular 4Z) spacetime, with a specified physical energy-momentum tensor, is an 
interesting and important question. However, it faces a fundamental problem. Namely, that the effective 4D equations 
for gravity contain a source term, that is the spacetime projection of the 5D Weyl tensor, which is unknown without 
specifying the properties of the metric in 5D [H]. Therefore, the best way of deriving some practical benefit from 
Campbell-Maagard's theorem is by solving the equations in 51?, under some general geometrical properties, and then 
study the possible interpretations as well as the characteristics of the matter sources in AD that can generate these 
properties. 

From a single solution in 5D one may construct a great variety of scenarios in AD ranging from static configurations 
to cosmological solutions [9], which reflects the fact that (a) there are many ways of embedding a AD spacetime in 
a given higher-dimensional manifold, and (b) that the effective equations is AD do not form a closed system of 
equations. Therefore, the study of exact solutions of 5D vacuum Einstein's equations is essential; not only because 
of the arguments (i)-(iii) given above, and for the formal mathematical aspects associated with the theory (e.g., the 
classification of 5D manifolds) but also, for the growing importance of the application of Kaluza-Klein to cosmological 
and astrophysical phenomena in AD. The study of all possible scenarios is important in order to be able to predict 
observable effects in AD caused by new physics from an extra dimension. 

Nearly all the solutions in 5D are obtained under the assumption of spatial spherical symmetry and the choice 
of one or more additional functional relations that simplify the field equations in such a way that they become fully 
integrable. The disadvantage of this approach is that different choices may lead to the same physical solution, in 
different disguise. 

In a recent work [10] we started a systematic investigation of the solutions of the 5D field equations Rab = for 
the case where the 5D metric possesses self-similar symmetry. We found the most general self-similar, homogeneous 
and isotropic, Ricci-flat cosmologies in 5D and showed that they can be interpreted, or used, as 5D Riemann-flat 
embeddings for spatially flat FRW cosmologies in AD. In this paper we continue that investigation, but we abandon 
the assumption of perfect spatial isotropy and consider cosmologies with spatial anisotropy. 

For the models to be physically realistic they have to closely resemble the standard FRW models. It is well known 
that, within the context of Bianchi classification, the generalization of the "flat", "open" and "closed" FRW models 
are cosmologies of Bianchi type I, V and IX, respectively [TT] . Cosmological observations today indicate that our 
universe is spatially flat |12|-[14j. Thus, Bianchi I models are the closest anisotropic approximation to the present 
stage of evolution of our universe. 

Consequently, in this paper we adopt the 5-dimensional line elemenlQ 

which on every hypersurface ip — constant, reduces to homogeneous and anisotropic Bianchi type-I cosmological 

^In this work we use the following conventions: t = x^, x = x^ , y = and z = x'^ are the usual spacetime coordinates; »/> = 
represents the coordinate along the "extra" dimension; the signature of the 5D metric is (+, — , — , — , e) where e can be either —1 or +1 
depending on whether the extra dimension is spacelike or timelike. The range of tensor indices is A, B... = — 4 and fi, u, ... = — 3. 
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models with flat spatial sections. A simplified version of this metric has been studied by Roque and Seller [15j 
under the assumption of separability. In this work we obtain the general solutions to the field equations Rab = 
under the assumption that the metric ([T]) possesses self-similar symmetry. This assumption is motivated not only 
by the fact that FRW models are self-similar in iD [TB] and 51? [TU] , but also by the "similarity hypothesis" of Carr 
and Coley [17 , which asserts that under a variety of physical circumstances, both homogeneous and inhomogeneous 
cosmological models (in 4D) naturally evolve to a self-similar fornjl. 

We demonstrate that the field equations in 5D allow three types of similarity solutions: homothetic, conformal 
and "wave-like" . We show how to integrate the field equations for the first two types of similarity. We find that the 
solutions (a) depend on one arbitrary function and three arbitrary parameters; (b) require the extra dimension to be 
spacelike, and (c) generalize the family of homogeneous and isotropic 5D models discussed in [10] . In accordance with 
Campbell-Maagard's theorem they can be used, or interpreted, as 5D embeddings for four-dimensional subspaces. 
Depending on how we do the actual embedding, one can construct AD anisotropic cosmological models that are either 
homogeneous or inhomogeneous, although homogeneous isotropic cosmologies are recovered as particular cases. 

By means of dimensional reduction of the metric in 5D, we analyze, in some detail, possible applications of the 
family of homothetic solutions in 5D. We find that all the models constructed in AD, exhibit some type of self- 
similarity. Those that are inhomogeneous and anisotropic inherit the homothetic symmetry from the 5D embedding. 
The rest of them are either partially homothetic or homothetic along some AD vector. Thus, all these models 
represent different self-similar scenarios in AD and, as such, can be relevant to the similarity hypothesis mentioned 
above. 

The paper is organized as follows. In section 2, we deduce the shape of the similarity variable. In section 3, we 
solve the field equations with homothetic and conformal symmetry in 5D. In section 4, we present some general 
properties of the solutions in 5D. In section 5 we show different admissible interpretations in AD, but an exhaustive 
treatment of all possibilities is beyond the scope of this work. In section 6 we give a summary. 

2 Similarity variable 

In a self-similar model, by a suitable transformation of coordinates all the dimensionless quantities can be put in a 
form where they are functions only of a single variable (say C) |22j-j29j. Thus, in the case under consideration, in 
"self-similar" coordinates {i,x,y,z,'4>), the line element ([1]) can be written as 

dS^ = e^'f^^dP - e^(^)dx2 - e^^^^Uf - e^^^^dz^ + ee^^^U^. (2) 
where C, is some function of t and ■)/', viz., 

C = C(i,^)- (3) 

The metric functions !^(C), A(C), m(C)i ^{C)i "^(C): well as C and the signature coefficient e, have to satisfy the field 
equations Rab = 0. In order to avoid misunderstanding, it is worth mentioning that in the literature the concept of 
self-similarity is frequently equated with homothetic symmetry. In this work we follow the traditional nomenclature 
used in [22]-[24]. Among other things, we will see that the self-similar line element ([2]) does not necessarily admit a 
homothetic Killing vector. In what follows we are going to suppress the bar over the self-similar coordinates. 

The shape of the similarity variable Q is determined by the field equations. Firstly, we note that from R^^ = it 
follows that 

2 (Ac + + ^c) C' + [(A^ + al) - (Ac + /ic + ac) [v^ + ^c) + 2 (Acc + A^CC + ^Cc)] C'C = 0, (4) 

where denotes derivative of / with respect to C,; dots and primes stand for derivatives with respect to t and tp, 
respectively. 

■^Although there are some counterexamples to this hypothesis |18) . |19j . there is a strong evidence that self-similar models play a 
significant role at asymptotic regimes I20| -|21 | . 
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2.1 Homothetic and conformal self-similar variable 



Equation ^ implies that, in the general case where C,' 0, the assumed self-similarity requires that the ratio 
[C'/(C'C)] be some function of (. Clearly, any separable function of t and ip will do the job. Therefore, without loss 
of generality we can set 

where T and F are some functions of their arguments. Consequently, the equations i?oo = 0; ^ii = Oj R22 — 0, R33 = 
and i?44 = they all have the structure 



f2\ 



fF\ 



-BH MiO+ — NiO+ 2- 



pi2 \ 



pi2 \ 



FF" 

pi2 



P{0 + QiC) = o, 



(6) 



where M , N, P and Q symbolize the corresponding functions of ^ in these equations. Therefore, in order to preserve 
the self-similar symmetry, we have to require 

FF" 

2-^='. (7) 
where I is a dimcnsionless separation constant. Integrating this expression we obtain 

F' = pF^ 

where p is a constant of integration. Thus, we find 





f2 




p2p{4-2l) 



C(^-2", and 



1} 
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m- 


p2p{3-2l} 



(8) 



(3-2/) 



(9) 



Consistency of ^ demands the quantities inside the square brackets to be constants, which requires 

tVC-i) for I ^ 1, 



T ' 



(10) 



for I = 1, 



where q is some constan10. Consequently, without loss of generality we can set 

I [^J 

c = 



(11) 



for / ^ 1, 
^) for I = 1. 

We will see that the variables ( for I ^ 1 and I = 1 lead to homothetic and conformal solutions in 51?, respectively. 



2.2 Wave-like self-similar variable 

Finally, in the particular case where C = 0, it is not difficult to show that the similarity variable must have the form 

C = c^oi + fcoV-, (12) 

where uq and fcp a-re some constants with the appropriate units. Thus, in this case the metric functions have a 
dependence of time and the extra coordinate like in traveling waves or pulses propagating along the fifth dimension. 
Planes waves or wave-like solutions in 5D, in the case of spatial spherical symmetry, have been studied by Wesson, 
Liu and Seahra [30], [31j; Horowitz, Low and Zee as well as by the present author [33] . 

^We note that p and q must have the units of (length)"^. 



4 



3 Field equations 



The field equations Rab — look quite different for C' 7^ and (' — 0. In what follows we will consider the general 
case where C,' ^ 0, for which the metric functions depend on pT|) . Wave-like solutions will be considered elsewhere. 

Let us notice that for I = 2 the similarity variable becomes C{i=2) = C = (^/V')- Thus, for an arbitrary Z ^ 1 we 
can write C = What this means is that we can first study the integration of the equations for I = 2, and 

then obtain the solution for any value of ^ by a simple transformation 



3.1 Integrating the field equations for / = 2 

The non-vanishing components of the Ricci tensor in 5D are 
1. Rqo — 0, 



2. i?ii = 0, 



3. R22 = 



4. -R33 = 



5. Raa = 0, 



2A, 



= 0. 



4ee'^ + (£e--^V)(u;5-;.5) 



= 0. 



e'e'') A^ + /i^ + (Tj + ^5 



2tJ: 



[4 {v^ + + + wj) + e (i^l + A| + + 4 + w|) 2^ {v^^ + A55 + /i« + + 
6. Finally, i?04 = 0, 

2^ (A« + + a«) + (2 - C^^e - (Aj + (75 + M?) + ^ (A^ + m| + ^^1) = 0. 

We note that equations p^ . and pT|) require 

A^t ^ ^ 

A? o-j " 

Therefore, without loss of generality one can set 

e^=Af"(0; e^'^Bff^iO, e'^^Cf^O, 



(13) 



(14) 
(15) 
(16) 
(17) 



(18) 



(19) 



(20) 



(21) 



where A, B, C are constants; / is some function of the variable ^ = (VV-"); a, [3 and 7 are arbitrary parameters. 
As a consequence, Ru — 0, R22 = and i?33 = reduce to 
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- i (j.^ - 0.5) (ee"-CV)+ 2^6^ = 0. (22) 



On the other hand Rq4 yields 

(« + /3 + 7) - - + 1^ + 2(^2 + + ^2 _ ^ _ ^ _ ^) 1^^^ ^ (23) 

from which we get 

where E is a, constant of integration, and 

a =(« + /? + 7), b={a^ + 13^ -a- (25) 

Clearly, from ((24)) it follows that there are different scenarios depending on whether [v + oj) or [v + oj) = 0. We 
now proceed to discuss them separately. 

3.1.1 General solution for / = 2 

For the general case where v + uj ^ in we substitute e"/^ = E^f^''/°'^f^e^'^/^ into ^ and Since these 

are second order differential equations for v, we obtain two differential equations containing /^jj, the third derivative 
of /(C)- Next, we isolate the /^^^ obtained from i?oo = and R44 = 0, say /^^^(-Roo) and /^^^(i?44), respectively. 
Then, from the compatibility condition 

MiRoo) = /«c(i?44), (26) 

we find 

g2^^_g^2^4^(2b/a)^2^ (27) 

This expression demands the extra dimension tp to be spacelike, i.e., e = — 1. In addition from ([24]) we find 



^^y(b/a)j^^ and e'^^) =C^e''(«). (28) 

Now, it is easy to verify that the field equation (|22p is identically satisfied for = £^"^6" and e = — 1. 
Collecting results, we have showed that the 5D line element 

= Ef^^'^^J^dt^ - Af^dx^ - Bf^dy^ - Cf^dz^ - E^'^ f^^'"^ f^d^^ , (29) 

is the general solution of the field equations (fT4l) -(fT9 |) . for any arbitrary function / — f{£) with £ = t/ip. 

3.1.2 Particular solution for / = 2 

For the case where j/ + = 0, but v = —uj ^ 0, integrating ([^^ we obtain the solution as follows 

dS^ = dt^ - Ah^^dx^ - Bh^l^dy^ - Ch^^dz^ - ^d^j^, with h=[C\\n£ + C2]"/(''+''), (30) 
where Ci and C2 are constants of integration. It should be mentioned that this is the unique family of solutions with 

iy=-UJ^O. 

In the case where v ^ —lo = the field equations Rab = require a = a + /3 + 7 = 0,i.e., h — 1. Thus, this case 
yields flat Minkowski space in 5D and AD. 
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3.2 Solutions for / ^ 1 

In order to construct the solutions with Z 7^ 1, we now use the transformation (jl3p . The second equation in 
becomes 



e 



'^(C) ^ (^2(/-i)g^(C)^ (31) 
On the other hand /e C^~'/c- Thus, from (|29p we obtain the general solution for Z 7^ 1 as 

dS^ = EC^-' f^dt^ - Af^'dx^ - Bff'dy' - Cf^dz^ - ^C'/^'/^Vc'^V''. (32) 
In addition, from ((30|l we get 

dS^ = (^) " - Bh^f^dy^ - Ch^^dz^ - C'^MV^ with h = [Ci InC + (72]"/^"+''', (33) 

which in terms of t and ?A is identical to pop . However, we put it here in explicit form in order to develop solutions 
with 1 = 1. 

3.3 Solutions for ^ = 1 

The field equations Rab — require p = q. Therefore the similarity variable now becomes 



In this case the solutions are readily obtained from ([5^ and ([55]) just by setting I = 1 and replacing ( by (. Thus, 
from PT|) we get 

e-(C) ^ e'f^^^ (35) 
Consequently, the general solution with Z = 1 is given by 

dS'^ = ECf^'^^fi^dt^ - Af^dx^ - Bff^dy^ - Cf^dz"^ - ECf^''/''^f^diP'^. (36) 
Likewise, from p3p we get the particular solution 

dS'^ = df - Ah^'^dx^ - Bh^l^dy^ - Ch'^^dz^ - dijj^ , with /i = [Ci InC + (72]"/^"+^'. (37) 

3.4 Solutions for a = a + /? + 7 = 

It should be noted that ^ and ^ reqmre a = (a + /? + 7) 7^ 0. If a 0, then from ^ it follows that 

either 6 = 0, or / = constant. In both cases the 5Z) metric becomes 

dS^ = e'^'^^^df - dx^ - dy^ - dz^ + ee'^^'^U^?. (38) 

At this point, it is crucial to emphasize that (I31|) and (I35p hold for any particular value of a,/? and 7. Therefore, 
the line elements 

/ \ i/C-i) 

dS^=e''^^Ue-dx^-dy^-dz^-e^'-^^e-^^U4,\ C = f ^ , (39) 

and 

dS^ = e'^^^Ue - dx^ - dy2 _ ^^2 „ eKC)rf^2^ ^~ ^ ' ^ ^^^^ 



constitute the most general solutions of the field equations Rab = 0, for the metric (|38p . We also note that (|30p . 
([55)) and ([57]) require (a + 6) = + + 7^ 7^ 0. This is always satisfied unless a = /3 = 7 = 0, in which case the 
line element has the form 



4 Some properties of the solutions 

The general solutions (j32p and p6p depend on one arbitrary function, of the corresponding similarity variable, and 
contain three arbitrary parameters a, (3 and 7. In addition, the field equations require the extra dimension to be 
spacelike (e = — 1). For a — f3 = j we recover cosmological models with spatial spherical symmetry. Indeed, setting 
= f^^ = P'^ = e'^'^^-'; changing coordinates: x — r sin 9 cos (j), y = rsmOsincj), z — rcosO; re-naming the 
constants, without loss of generality (|32[) can be written as 



dS^ = (^i^ C^-^-'hce^^<'>^' - e^K) [dr^ + {d9^ + sin' Od^'')] - C'Ace^^^^/^ d^'- (41) 



Likewise (l36l) becomes 



These two metrics represent the most general self-similar, homogeneous and isotropic, Ricci-flat cosmologies in 5D 

m- 

Although solutions (j36|) are formally obtained from ((32|) just by setting I — 1, they have different geometrical 
properties. In particular, solutions (I32p admit a homothetic Killing vector in 51? for any values of a, (3 and 7, namely, 

C^g^P = 2glf\ with Cfi^^^ = {t, x, y, z, V), (43) 

where g^^^"* is the metric ([5^ and denotes the Lie derivative along the 51? vector Cfi^i)- the other hand, 
solutions (|36p are self-similar but do not admit a homothetic Killing vector, except for a = /3 = 7. In general, they 
admit an infinitesimal conformal transformation parameterized by some function H = H{t + tp), viz., 

Ccg'i';'^ =2Hg'i='\ with C(?=i) = (H, xH, yH, zH, H), and H = H{t + ij). (44) 

Another difference is that (|32|1 is invariant under the transformation: t iijj, it (double Wick rotations) 

and E —E. Indeed, it is easy to verify that 

dS"" = EC}f"''^f^de - Af"dx' - Bff'dy' - Cf^dz' - EC^'-'^ f^'/^^ f^d^p' , with C = ^ (45) 

is also a solution of the field equations. The same occurs with but the metric functions in the "Wick rotated" 
line element are now complex functions. 

However, the solutions share some important properties. For example, one can show that there are only two 
cases where the components of the 51? Riemann tensor vanish: (i) / = constant (or a = /3 = 7 = 0), and (ii) 
a — P — ^ 0. In the first case the line element and the metric functions are given by ((38|) . In the second 
case the spatial sections t ~ constant and ^ = constant are flat and possess spherical symmetry 10]. In any other 
circumstance, the manifold is curved in 51?. As an illustration we present here the Riemann tensor in 51?, calculated 
with the line element (l29l). 
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-R0II4 = ^^0101 = 
-R0224 = C^0202 — 

-R0334 = ^^0303 = = — 75 -■ (46) 

Finally, it is essential to mention that the solutions discussed in the preceding section are invariant under the change 
(x, y, z) ^ ip, which is a consequence of the fact that the extra coordinate '0 is spacelike. For example, if we change 
■0 <-> z and denote 

F = F{r)), with ?7 = -, (47) 
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then the hne element 

dS^ = EF^^/'^'^F,^ [dt^ - ifdz^) ~ AF'^°'dx^ - BF^'^dy^ ~ CF^^diP^, (48) 

which we have constructed from (|29p. is also a solution of the field equations Rab — 0. Certainly, the same is true for 
all solutions discussed in section 3. It is important to note that the metrics obtained under the change (x, y, z) ip 
(i) are independent on the "extra" dimension -0, and (ii) the spatial sections defined by t = constant, ip = constant 
are non-flat. In the next section we will see that this symmetry is crucial for the interpretation of the solutions in 
4D. 



5 Interpretation in AD 

We have already mentioned that there are different ways of producing, or embedding, a AD spacetime in a given 
higher-dimensional manifold. However, the most popular approach is based on three different assumptions. First, 
that we can use the coordinate frame [51]. Second, that our AD spacetime can be recovered by going onto a 
hypersurface : ^ = ipo — constant, which is orthogonal to the 5D unit vector 

n^ = 4^, UAU^^e, (49) 



along the extra dimension. Third, that the physical metric of the spacetime can be identified with the one induced 
on E^. 

For a line element of the form 

dS^ ^ 9^^{xP, ij)dx''dx'' + e<i>^{xP, ip)dip^ , (50) 

the induced metric on hypersurfaces E^ is just g^^^ i.e., the AD part of the metric in bD. The crucial moment 
is that, although the energy-momentum tensor (EMT) in 51? is zero, to an observer confined to making physical 
measurements in our ordinary spacetime, and not aware of the extra dimension, the spacetime is not empty but 
contains (effective) matter whose EMT, '-^-'Tq,^, is determined by the Einstein equations in AD, namely 

= 8^ (4)T„^ = -e {K^y^K^p - K^K^p) + ^g^p {KxpK^" - [K^f) - eE^p, (51) 

where K^^ is the extrinsic curvature 

K^p = -Uga0 = (52) 
Efj^i, is the projection of the bulk Weyl tensor '^^^Cabcd orthogonal to n"^, i.e., "parallel" to spacetime, viz., 

1:1 (5) /-I -A-B 1 dKafi p 'I'q;/? ,.„^ 



and $Q = 9$ /dx" . It is important to mention that the effective matter content of the spacetime is the same whether 
we interpret it in space-time-matter theory j35|-[36j. or in a Z2 symmetric brane universe |37j . Indeed, these two 
theories are mathematically equivalent, although they have different motivation and physical interpretation [38| . 



5.1 Homogeneous cosmological models 

Since C = it follows that the models with I ^ 1, discussed in Section 3.1, are physically indistinguishable 

from those with I ~ 2. The only reason for considering them here is a mathematical one: they allow us to generate 
non-homothetic 5D models (^ = 1) from homothetic ones (^ 7^ 1). In what remains we focus our attention to metrics 
and which admit a homothetic Killing vector in 51?, namely and therefore are relevant to the similarity 
hypothesis mentioned in the Introductior0. 

''We shall come back to discuss the 4-dimonsional interpretation of the non-homothetic 5D models (/ = 1) elsewhere. 
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However, by virtue of the symmetry of the solutions in 5D, which in the case under consideration is expressed by 
and P5)l . there are a number of distinct possible scenarios in AD. An exhaustive treatment of all these scenarios 
is beyond the purview of this work. Rather, in this section we restrict ourselves to presenting some representative 
models. 

In the approach under consideration, the spacetime metric generated by (j29p is given by 

ds^ = dSf^^ = gf.^dxf'dx'' = E f'^^^''^ f^dt^ - Af^'dx^ - Bf^dy^ - Cf'^dz'^. (54) 

Let us emphasize that, on each hypersurface S^, this metric is a function of time only. Therefore, ([5^ yields 
homogeneous, anisotropic cosmological models of Bianchi type 1. 

A general feature of metrics ([5^ is that they are partially honiothetic, i.e., 

^cgoo = 0, Ct^gij = 2gij, with = (0, x, y, 2), (55) 

for any arbitrary function /. 

5.1.1 Models with homothetic symmetry in AD 

The condition of homothetic symmetry in AD singles out a unique family of metrics, which are generated by the 
choice / = ^'^/ ("■+''') (See the Appendix). The homothetic AD vector as well as the embedding solution in 5D are 
given by (|A-4[) and (jA-5[) . respectively. The AD line element with homothetic symmetry is given by (jA-6[) . 

a + /3 + 7 
a + b a2 + /?2 + ^2 

We note that the volume of the spatial slices goes like t'^'^. Since afc > 0, it follows that the volume of space is 
increasing from an original big-bang at t = 0. This is regardless of whether some particular direction is contracting. 
The effective EMT ^ 

g^^o^fc^(a^+_J7+M, Ti = n.r°, Ti = nyT^, T^^n^T^, (57) 

where nx, riy and Uz are constants given by 

_ {-a + f3 + j) _ {-(3 + g + 7) _ (-7 + a + /3) 

~ a + /3 + 7 ' "'^ ~ a + /5 + 7 ' " a + /5 + 7 
We note that 

Tlx + Uy + Uz = 1, (59) 

for any value of a, (3 and 7. Let us notice some particular cases: (I) If one of the coefficients vanishes, say a = 0, 
then Ux = ^ and Uy = — n^; (II) If two coefficients are equal, say a = /3, then Ux ^ ny. In particular, if a = P ~ 0, 
then Ux = ny — I and Uz = —1; (III) If a = /3 = 7 ^ 0, then Ux — Uy — Uz — 1/3. The homothetic model exhibits a 
number of interesting properties, however a detailed study of the matter sources that satisfy ([57)) -(|5 ^ is outside the 
scope of the present work. 

5.1.2 Perfect fluid with stiff equation of state 

The function /(^) in (|54p can be determined by imposing some "physical" conditions on the effective picture in AD. 
If we assume that the effective matter behaves like a perfect fluid, then we obtain 

a + /3 + 7 
+ /32 + 72 + 4(a/3 + ^7 + /37) ' 

and Ci, C2 are some constants. The effective density p ~ and pressure p = —T^ = —T2 ~ — T| in AD satisfy the 
stiff equation of state, namely, 

_ [7(a + + ap] m Ci 



(58) 



/(0 = (Cie + C2)'", where , , 2 T . , , (60) 



^To simplify the notation, in what follows we suppress the index in '^'T, 



a/3- 
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5.1.3 Homothetic 5D embedding for the 4D Kasner spacetime 



For completeness, we should mention that on each spacetime section we have 

Snn,^l^^^±fl±^^li)\ (62) 
Therefore, vacuum solutions in AD require 7 ~ ^af3/{a + (3). Then, from the vanishing of we obtain 

/(0-(Cie + C2)°^+i+''^ , (63) 

where Ci and C2 are constants of integration. Now, setting Ci — [{a^ + af3 + 0^)/E{a + /?)], we obtain 500 = 1- 
Thus, the metric in bD becomes 



dS^ = - A{Cii + C2fP'dx^ - B{Ci^ + C-ifP^dy^ - C{Ci^ + C2fP'dz^ - (^^^dil)'^, (64) 

where 

■ P4 = 1- (65) 



a{a + P) Pia + (3) a(3 + a(3 + (3^ 

We note that 

Pi+P2+P3^i, pI+pI+pI^^- (66) 

Therefore, can be interpreted as a homothetic 5D embedding for the 4D Kasner spacetime. It is important to 
point out two things. Firstly, that although (j64p represents an empty space in 5Z?, it is not equivalent to the usual 
5-dimensional Kasner spaceO Secondly, it illustrates that the effective EMT in 4Z3, which is given by (|5ip can vanish 
even though the induced metric in AD explicitly depends on the extra coordinate. This is totally different from the 
isotropic cosmological models, for which the effective EMT is always nonzero when the 5D metric has a explicit 
dependence on the extra coordinate. In fact, this is a general attribute of models with spatial spherical symmetry 



5.2 Inhomogeneous models: nonstatic vacuum solutions 

Let us now consider the 5D line element ((48]). In the approach under consideration, the spacetime metric is given by 
ds^ EE dS^^^ = g^^dx^'dx" = EF^^'^^Frj [dt^ ~ ifdz^) ~ AF^'^dx^ ~ BF^f^dy^, (67) 

where is a function of the self-similar variable 77 = t/z. These models present two main characteristics. First, they 
show homothetic symmetry along the AD vector 

= (i, X, y, z), (68) 

where the subindex p indicates that is the spacetime projection of the 5D homothetic vector defined in (|43p . 
Second, since the spacetime part of the generating 5D metric (|48p is independent of ip, it follows that K^j^^, — and 
consequently the effective EMT is 

SttT^. = S^. = ^|^, (69) 
we recall that the general solution requires e = — 1. Therefore, 

T = To° + + T| + = 0, (70) 
which is a consequence of the fact that E'pj, is traceless. 

^We recall the reader that a Kasner space in 5D is given by dS^ = dt^ - At^Pirfx^ - Bi^P^dy'^ - Ct^P^dz'^ + tEfiPidip^ with 
Pi +P2 +P3 +P4 = 1, +P2 +P3 +P4 = 1- 
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5.2.1 Solutions in conventional AD general relativity 

The nonvanishing components of the EMT are 







(f) 



7 [7 (a + /3) - - /?■ 
87r(a + P + -f) E 



'21 r 



^2 ^(2a+6)/a 



(71) 



They can be interpreted as "pancake-hke" distributions of matter, with energy flow along the direction of symmetry 
z. A detailed study of the matter distribution ((7T|) is beyond the scope of the present work. 
For another interpretation, we note that T^^ = when 



a + P 



a ^ -p. 



(72) 



With this choice, the line element (j67p yields a two-parameter family of non-static inhomogeneous vacuum solutions 
in conventional AD general relativity. For these solutions, the non-vanishing components of the Riemann tensor in 
5D and AD are all proportional to (a — f3)Frj. Consequently, for a ^ /3 the vacuum solutions given by ((67|) and ((72|) 
are not equivalent to a AD Minkowski spacetime and cannot be embedded in a Riemann- flat manifold in 5D. 



5.2.2 Vacuum solutions in the Randall-Sundrum (RS2) braneworld scenario 

Without going into deep technical details, in the Randall & Sundrum braneworld scenario |39| our universe is 
identified with a fixed singular hypersurface (the brane) embedded in a 5-dimensional bulk with Z2 symmetry 
with respect to S^^[37]. Due to the presence of matter on the brane, which is described by an EMT that we denote 
as r^^, the extrinsic curvature K^i, is discontinuous across Israel's boundary conditions |41J relate the jump of 
K^^ to r^j,, namely, (i^^,y|j^+ — if^^ij.- ) — A:(t^^ — hTg^iv): where A: is a constant with the appropriate units. Next, 

Z2 symmetry implies K^^,^+ — -~K^^^^~ = K^i, . Therefore, when the spacetime (the brane) is empty (r^^ = 0) 

'*b 

it follows that K^^^ ~ 0. Thus, in braneworld theory the vacuum field equations are obtained from ([5T|) as 

(4)G„^ = Stt (4)t„;5 = -eE^p, (73) 

which means that, in braneworld theory, the effective geometrical matter is AD is traceless. Therefore, any solution 
of general relativity with T = is a solution of the vacuum equations in braneworld. 

The conclusion is that ([67)1 can be interpreted as vacuum solutions in the Randall-Sundrum braneworld scenario. 



6 Summary 

Self-similar (homothetic) symmetry seems to be very important in naturt0. In cosmological applications, within the 
context of conventional general relativity in AD, there is a strong evidence that many homogeneous and inhomoge- 
neous cosmological models can be approximated by self-similar models in the asymptotic regimes, i.e., near the initial 
cosmological singularity and at late times |17j-[21j. In order to avoid misunderstanding, it is important to emphasize 
that in the literature the concept of self-similarity is frequently equated with homothetic symmetry. In this paper we 
have followed the traditional terminology used in [22]-[2l|, where self-similarity means that the field equations are 
functions only of a single variable, which in turn allows to reduce them to a system of ordinary differential equations. 

Observations indicate that on large scales (^ 100 Mpc) the universe is homogeneous and isotropic and well 
described by spatially-flat FRW cosmologies, which are self-similar in AD and hD. However, there is no reason 
to expect such features at the early stages of the evolution of the universe. Rather, it is generally accepted that 

^We would like to share with our readers a beautiful quotation about self-similar (homothetic) symmetry, ascribed to Manfred 
Schroeder, that we learned from a paper by Adrian Popesku |42| : "The unifying concept underlying fractals, chaos and power laws is 
self-similarity. Self-similarity, or invariance against changes in scale or size, is an attribute of many laws of nature and innumerable 
phenomena in the world around us. Self-similarity is, in fact, one of the decisive symmetries that shapes our universe and our efforts to 
comprehend it." 
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anisotropy could have played a significant role in the early universe and that it has been fading away in the course 
of cosmic evolution. Therefore, the study of cosmological models that are anisotropic and self-similar appears to be 
of primary interest. 

In this work we have investigated such cosmological models within the context of theories of Kaluza-Klein type, 
with a large extra dimension. The ladder to go between the physics in AD and 51? is provided by Campbell-Maagard's 
embedding theorems, which guarantee that any solution of the AD Einstein equations of general relativity may be 
embedded in a solution of the 5D vacuum Einstein equations. 

For the case where the 5D metric is diagonal and is a function only of time and the "extra" coordinate ^p, we have 
shown that there are three possible forms for the similarity variable: £, = t/ip; ^ = {e* /e^Y^ and C, = ivot + kotp. They 
correspond to three different physical situations, which are described by homothetic, conformal and plane wave-like 
solutions in 5D, respectively. Thus, in the traditional terminology the concept of self-similarity includes conformal 
and wave-like solutions. 

In section 3, we have constructed the most general homothetic and conformal solutions in 5D to the field equations 
(|14p - (jl9p . They are given in terms of one arbitrary function, of the appropriate similarity variable, and three arbitrary 
parameters, viz., a, (3, and 7. However, the situation is completely different for the case where the bD metric is wave- 
like. This case is more involved and requires a separate discussion. 

Following Campbell-Maagard's embedding theorems, the connection to AD is deduced after choosing an embed- 
ding. We have found that the solutions that are homothetic in bD are relevant to the similarity hypothesis because 
all the models constructed in 4Z?, by means of dimensional reduction of the metric in 5-D, exhibit some type of 
self-similarity. Specifically, we have seen that those that are inhomogcneous and anisotropic inherit the homothetic 
symmetry from the 51? embedding. The rest of them are either partially homothetic or homothetic along some AD 
vector. 

Our discussion in section 5 illustrates two important things. Firstly, that when the extra dimension is spacelike 
the concepts of AD spatial homogeneity and AD spatial flatness become totally dependent on the embedding, i.e., 
on how the coordinates for our 3-dimensional space are chosen. Secondly, that the interpretation in AD crucially 
depends on the theory under consideration. In consequence, much work is still needed in order to understand AD 
physical models as Lorentzian hypersurfaces in pseudo-Riemannian 51? spaces. 

In summary: our work (I) generalizes the 5D embeddings used for the FLRW models; (II) shows that anisotropic 
cosmologies are, in general, curved in 5D, in contrast with FLRW models which can always be embedded in a 5D 
Riemann-flat (Minkowski) manifold; (III) reveals that anisotropic cosmologies can be curved and devoid of matter, 
both in hD and AD, even when the metric in 5D explicitly depends on the extra coordinate, which is quite different 
from isotropic cosmological models where the effective EMT is always nonzero when the 5D metric has a explicit 
dependence on the extra coordinate. In fact, this is a general attribute of models with spatial spherical symmetry 

m- 

To finish this paper, we should mention that from a mathematical point of view the shape of our solutions can 
be simplified if we introduce "null- like" coordinates (We thank Philippe Spindel for pointing this out 03]). As an 
illustration, let us consider ((29l) . The {t — ip) part of that metric can be factorized and written as E f^/°- f(^{dt ~\- 
idil;){dt - idip) = 2dudv. Then, setting u = f{°-+b)/o- the hue element becomes 

dS"^ = 2dudv - AuP'dx^ - Bu^^dy^ - Cu^'dz^, (74) 

where pi = 2aa/ {a + b), p2 — 2/3a/{a + b) and ps — 27a/ (a + b). We note that {pi — 1)^ — 3. Although in this 

coordinates the solution looks mathematically simpler than in the original form, its physical interpretation in AD is 
elusive since it not clear how to define the hypersurfaces orthogonal to the 5D unit vector along the extra 
dimension. Therefore, in these coordinates neither the dimensional reduction discussed in Section 5 nor Campbell's 
theorem can be used for the AD interpretation of the solutions. 

This investigation can be extended, or generalized, in different ways. They follow from the fact that we have not 
fully examined the possible AD interpretations of the 5-D-homothetic solutions ([5^ . Neither the conformal solutions 

^As it was mentioned in the Introduction, the physicaUy interesting case of finding a 5D embedding for a given AD spacetime, with a 
specified physical EMT, is a very difficult task because the effective 4D equations for gravity contain a source term, that is the spacetime 
projection of the 5D Weyl tensor, which is unknown without specifying the properties of the metric in 5D 8 . 
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(|36p . nor the plane wave-like solutions induced by the similarity variable (jl2p . have been studied. These are important 
topics that should be addressed. 



Appendix: Homothetic symmetry on 

Our aim here is to show that the requirement of homothetic symmetry on singles out one specific metric in AD. 
With this aim, let us take the Lie derivative of ((54)) along the AD vector 

= (Co t, Ci X, C2 y, C3 z) , (A-1) 

where h stands for "homothetic" and Co, C2, C3, C3 are some constants. We obtain, 



3i ^3 

2goo Co 
2511 

2ff22 
2.g33 



1 + + 
2 U/ /« 



Ci + aC„C [ ^ 
C2+/3Coe( ^ 



C3+7Co^ ( ^ 



The condition Lq^q^i, = 2g^i, requires 

1 - Ci 1 - C2 1 - C3 



fiO-^'' and Co = 



2a 



aCo /3Co 7C0 "'"^ ^" a + k{a + b)' 

where k is some constant. Without loss of generality we can set Cq = 1, which implies 



k = 



-, and — [t, (1 ~ ka)x, (1 — k(3)y, (1 — kj)z] 



Thus, the 5D metric 



generates, on every E^, the AD line element 



ds^ EE dS'Pv 



- A e'^^dx^ - B e^'^dy' - c e^'^dz', 



2/3fc^„2 



(A-2) 
(A-3) 

(A-4) 
(A-5) 
(A-6) 



which shows homothetic symmetry along Qf^. Here we have introduced the constants A, B and C for dimensional 
consistency. It should be noted that C^f^ is not the spacetime projection of the 5Z) vector C(^^i) — {t, x, y, z, tp) 
defined in (03]). 
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